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ABSTRACT 

We have found that the two-point correlation function of the APM clusters of galaxies has an amplitude much 

• higher than that claimed by the APM group. As the richness limit increases from 1Z = 53 to 80, the correlation 

CNI \ length increases from 17.5 to 28.9 /i~ 1 Mpc. This indicates that the richness dependence of the APM cluster 

correlation function is also much stronger than what the APM group has reported. The richness dependence can 

J> ■ be represented by a fitting formula r$ — 0.53d c + 0.01, which is consistent with the Bahcall's formula ro — 0Ad c . 

CNI \ We have tried to find the possible reason for discrepancies. However, our estimates for the APM cluster correlation 

F"tT function are found to be robust against variation of the method of calculation and of sample definition. 
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■ I. INTRODUCTION 

Rich clusters of galaxies have been used for studies of large-scale structure on scales ~ 10 — 100 /i _1 Mpc. The 



spatial two-point correlation function for the Abell clusters has been estimated by many authors (see the review by 
Bahcall 1988; also Huchra et al. 1990; Posti: 
consistent in shape with the power law form 



^ | Bahcall 1988; also Huchra et al. 1990; Postman et al. 1992; Peacock & West 1992), which has been found to be 



U(r) = 7 (1) 

with the correlation length ro ~ 20 — 26 7i _1 Mpc and with the power law index 7 ~ 2. More recently new catalogs 
of clusters have been obtained by automated selecting algorithm from the Edinburgh-Durham Southern Galaxy 
Catalog (Lumsden et al. 1992), and from the APM Galaxy Survey (Dalton et al. 1997). Redshift surveys of these 
new clusters (Collins et al. 1995; Dalton et al. 1994a) have been used to estimate the spatial two-point correlation 
function (Nichol et al. 1992; Dalton et al. 1994&). The correlation length ro measured for these new clusters is 
reported to be lower than that for the Abell clusters, that is 14 < ro < 16 /i -1 Mpc. 

It has been argued that the counting radius (r c = 1.5 /i _1 Mpc) for the Abell clusters is so large that the 
catalog contains projection effects, which cause artificial line-of-sight correlations (Sutherland 1988; Efstathiou et 
al. 1992). In addition to this, it has been recognized that the intrinsically subjective nature of the Abell catalog 
causes problems in homogeneity and statistical completeness (Postman et al. 1986). However, Bahcall & West 
(1992) have suggested that the discrepancy between the Abell clusters and the APM clusters can be accounted by 
the richness dependence of cluster correlation amplitudes (Bahcall 1988; Bahcall & West 1992) 



r = 0.4d c , (2) 

where d c is the mean intercluster separation (Here d c — n c , where n c is the mean space density of clusters). 
But Croft et al. (1997) have analyzed richness subsamples of the APM clusters and argued that there is only a 
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weak dependence of correlation amplitude with cluster richness, and that this disagrees with results from the Abell 
clusters. 

We describe our statistical richness subsamples drawn from the APM cluster catalog in Section II. The method 
of calculation of the spatial two-point correlation function and our results are presented in Section III. In Section 
IV we discuss many possibilities which might affect the estimate of the correlation function. 



II. THE APM CLUSTER CATALOG 

Dalton et al.(1997) have published the catalog of APM clusters which is complete over richness range TZ > 40 
and the characteristic magnitude range 17.5 < m x < 19.4 (the distance estimation m x corresponds to Abell's 
mio )• The catalog covers the region of the sky 2l h < a < 5 h and -72?5 < S < -17?5 of the APM Galaxy 
Survey (Maddox et al. 1990; see also Loveday et al. 1996) and contains 957 clusters (see Fig. 1). The exact APM 



survey field definition can be found in the Astronomical Data Center ( http://adc.gsfc.nasa.gov/ ). Fig. 1 shows the 
distribution of 957 APM clusters (open circles), 185 APM fields (open squares), and 1456 holes (filled squares) in 
the (a, S) plane (Loveday et al. 1996). Holes are the regions contaminated by big bright objects. 

Although the 957 sample is complete in itself, it contains only 374 clusters with measured redshift. In the 
range of TZ > 53 and 17.5 < m x < 19.2, however, the completeness of the redshift sample is 90.9% and increases up 
to 100% at higher richness limits. So we generate four statistical subsamples with TZ > 53, TZ > 60, 1Z > 70, and 
1Z > 80, which contains 213, 149, 76, and 43 clusters, respectively. We have further restricted our APM samples to 
the declination range —65° < <5 < —25° to reduce the strong dependence of number of clusters in declination (see 
below). To convert the redshifts to comoving distances we assumed the Einstcin-de Sitter universe with fig = 1. 



III. CLUSTER CORRELATION FUNCTION 

We estimate the spatial two-point correlation functions using the Hamilton's estimator (1993) 

(DD)(RR) AN c N r 
tcclrJ " (DRY {N c -l)(N r -1) ' [ > 

which is evaluated to be less affected by uncertainties in the selection function for £ cc < 1 (Dalton et al. 19946). In 
equation (3), DD is the number of pairs in the sample with iV c clusters, RR is the number of pairs in a random sample 
with N r random points, and DR is the cluster-random pair count. The numerical factor 4 in the normalization 
term of pair counts accounts for the fact that we count each DD or RR pair only once. The uncertainties of £ cc are 
computed from the simple equation (5£ cc = (1 + £ CC )/V DD, which is easy to calculate but may underestimate the 
cosmic variance in comparison to simulations (Croft & Efstathiou 1994). 

Croft et al. (1997) have estimated the correlation function of the APM clusters located in the distance range of 
50 /i _1 Mpc < r < 500 h~ 1 Mpc. They used a selection function obtained by smoothing the distribution of clusters 
in redshift space over 40/V2 /i -1 Mpc by a Gaussian filter (dashed lines in Fig. 2). The random catalog is generated 
in accordance with this radial selection function within the sample boundaries. Equal weight is given to clusters 
and random points in the calculation of correlation functions. 

We have calculated the correlation function of the APM clusters in the exactly same way that Croft et 
al. (1997) have followed. Our estimates of correlation functions are plotted in Fig. 3 for four subsamples; the 
TZ > 53 (filled dots), TZ > 60 (triangles), TZ > 70 (squares), and TZ > 80 (open circles) samples. The cor- 
relation lengths are 17. 5l^g, 23.4^2'q, 25.1*37, and 28.917'j h~ 1 Mpc, respectively. For a comparison the fit- 
ting lines to Croft et al. (1997) 's TZ > 50 and TZ > 70 samples are drawn as solid and dashed lines, re- 
spectively. It can be seen that there are clear discrepancies between our results and Croft et al (1997)'s. 
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Fig. 1. The distribution of 957 APM clusters (open circles), 185 APM fields (open squares), and 1456 holes (filled squares). The long 
dashed lines represent the galactic latitude. 
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Fig. 2. Distributions of clusters with distance for subsamples with 1Z > 53 and TZ > 70. Short dashed curves are distributions 
smoothed by a 40/\/2 h~ 1 Mpc Gaussian. Solid lines are the selection function obtained by eq. (4). 



Amplitudes of correlation functions calculated by us are much higher both for 1Z > 53 and 1Z > 80 samples. 
Furthermore, the richness dependence is much stronger. The richness dependence of correlation functions of our 
APM subsamples can be described by a formula tq = 0.53d c + 0.01. This richness dependence is stronger even than 
that reported for the Abell clusters (Bahcall & West 1992); rg = 0Ad c , quite contrary to Croft et al. (1997)'s claim. 



IV. DISCUSSION 



We have looked for various possibilities that could explain the discrepancies between our and Croft et al. 
(1997) 's estimates of the APM cluster correlation function. 
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Fig. 3. Correlation functions of the APM cluster subsamples with 1Z > 53 (filled circles), TZ > 60 (triangles), TZ > 70 (squares), and 
1Z > 80 (open circles). We also present two fitted lines of Croft et al. (1997) for comparison; the solid line is for their sample B with 
rZ > 50, and the dashed line is for their 1Z > 70 sample. 



First, we changed the method calculating the selection function. Instead of using the distribution of clusters in 
redshift space smoothed over 40/v / 2 h~ 1 Mpc, we calculated the selection function from the conventional formula 



where is the solid angle of the SlirVGy region. &max,i cincl Drnin.i 

are the maximum and minimum distances 
to which the i-th cluster could be included in the sample for the faint and bright magnitude limits of the survey 
miim,upp = 19.2 and mu m ,low = 17.5, respectively. The selection functions calculated by this formula are shown in 
Fig. 2 for the 1Z > 53 and 1Z > 70 samples. We have calculated the correlation functions of the APM subsamples 
using these new selection functions. However, there is practically no change in the amplitude and in the richness 
dependence of the cluster correlation functions. 

Second, we limited the distance range of clusters to 170 /i _1 Mpc < r < 330 /i _1 Mpc instead of 50 /i _1 Mpc < 
r < 500 ft _1 Mpc to eliminate the possible dominance of shot noises. Again this did not affect our results. 

Third, we changed the way to give weights to clusters and random points in the calculation of correlation 
functions. Instead of giving equal weights we gave weights equal to the inverse of the selection function. To reduce 
the shot noise we limited the subsamples to the distance range of 170 h~ 1 Mpc < r < 330 /i _1 Mpc. This time, 
amplitudes of the correlation functions became slightly higher than those shown in Fig. 3. But this scheme still 
gives results which are inconsistent with Croft et al. (1997) 's weak cluster correlation function. 

Finally, we inspected the distribution of clusters in declination and galactic latitude spaces. If the sam- 
ple is affected by large galactic obscuration or airmass variation across the sky, a selection function varying 
on the sky should be taken into account. Fig. 4 shows the number of clusters in declination (squares) or 
in galactic latitude (circles) strips which have equal areas on the sky. Points are located at the centers of 
strips. It can be shown that, while there is no large fluctuation in the number of clusters in the galactic lati- 
tude space, there is a rather strong variation in the declination space for the 1Z > 53 and 1Z > 70 subsamples. 
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Fig. 4. Variation of the number of clusters with declination (squares) and galactic latitude (circles). 

Large number of clusters are concentrated in the declination strip centered at 8 = —48°. We first varied the decli- 
nation limits from the original range — 72?5 <J S <; — 17?5 to narrower ranges in the subsample with TZ > 53. The 
correlation function at short separations did not change but became somewhat steeper. However, the decrease in 
the amplitude of the correlation function at large separations is small, leaving the correlation function still within 
the la error limits. We then used the smoothed distribution of clusters in declination to generate the random 
catalog. This eliminates the effect of the strong variation of number of clusters with declination. The resulting 
correlation function has again somewhat lower amplitude at large separations, but the changes are again within the 
error limits. 

All these tests indicate that our estimate for the APM cluster correlation function is robust. It remains to be 
seen what has made APM group's results different from ours. 
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